Effect of a marginal inclination on pattern formation 
in a binary liquid mixture under thermal stress. 
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Convective motions in a fluid layer are affected by its orientation with respect to the gravitational 
field. We investigate the long-term stability of a thermally-stressed layer of a binary liquid mixture 
and show that pattern formation is strongly affected by marginal inclinations as small as a few 
milliradians. At small Rayleigh numbers the mass transfer is dominated by the induced large scale 
shear flow, while at larger Rayleigh numbers it is dominated by solutal convection. At the transition, 
the balance between the solutal and shear flows gives rise to drifting columnar flows moving in 
opposite directions along parallel lanes in a Super-Highway configuration. Our findings suggest a 
potential relevance for the experimental modeling of thermohaline convection in the oceans. 



PACS numbers: 47.20.Bp, 47.54.-r, 47.54.-r, 92.10.af 

Pattern formation in non-equilibrium systems arises 
from symmetry breaking of an isotropic initial state [T]. 
Whenever two symmetry-breaking mechanisms coexist 
their competition gives rise to a rich phase diagram. A 
typical example is represented by an inclined layer of liq- 
uid under the action of a temperature difference. Tilting 
the layer can determine a large scale shear flow (LSF). 
For fairly large inclinations theoretical and experi- 
mental studies [SHH] reveal rich spatio-temporal dynam- 
ics in the phase diagram, characterized by the presence 
of Busse oscillations, subharmonic oscillations as well as 
longitudinal and cross rolls. The transitions between dif- 
ferent regimes occur at angles of the order of tens of 
degrees, and pattern formation does not appear to be 
influenced by small inclination angles of the order of a 
few degrees or smaller. By adding a second component 
to the mixture, the compositional stratification also con- 
tributes to the convective behavior. A remarkable exam- 
ple is represented by thermohaline circulation in oceans, 
where both the local salinity and temperature of water 
contribute to the convective motions, the thickness of the 
layer of water being modulated by the seafloor [lOj [11] . 
The thermal stress applied to a liquid mixture via a tem- 
perature gradient can be quantified by the dimensionless 
Rayleigh number Ra — gApxh 3 /(j]Dt) [H], where g is 
the gravity acceleration, ApT the density difference gen- 
erated by thermal dilation of the liquid, h the liquid layer 
thickness, r\ its shear viscosity and Dt the thermal dif- 
fusivity. The presence of the thermal stress determines 
a separating flux of the two components and, in turn, 
a solutal density difference Ap c throughout the sample, 
whose value and orientation is quantified by the Soret 
coefficient St [Ml [14] . In the case of a mixture with pos- 
itive St, as the one used in our study, both ApT and Ap c 
contribute to destabilize a sample heated from below. 

In this letter we investigate the influence of a small 



inclination angle in the range mrad< a < 36 mrad 
(a < 2°) on the long-term stability of a layer of a binary 
liquid mixture under thermal stress below the thresh- 
old for Rayleigh-Benard convection (Ra c ~ 1700). This 
model system mimicks many of the features of thermo- 
haline circulation. Interestingly, we find the stability of 
the mixture to be strongly influenced already by small 
inclination angles, in contrast to the behaviour of simple 
liquids. The observed convective patterns are arranged 
into a peculiar Super-Highway configuration, reminiscent 
of the traffic of vehicles during rush hours. Such pattern 
differs dramatically from the square, roll, and cross-roll 
patterns previously reported experimentally [151 116| and 
theoretically [T71 [T5] for binary liquid mixtures heated 
from below in the absence of inclination. 

The binary mixture selected for this study is 
isobutylbenzene/n-dodecane at 50% weight fraction. The 
choice of this mixture is motivated by the availability of 
an accurate and extensive set of thermophysical proper- 
ties [TS]. The sample is delimited by two 8 x 40 x 40 
mm 3 square sapphire windows kept at fixed distance and 
sealed by a circular Viton O-ring gasket with an inner 
radius R = 13 mm. The geometry corresponds to a 
moderately high aspect ratio R/h = 10 of the sample. 
The sapphire windows are in contact with two annu- 
lar thermo electric devices connected to two independent 
proportional-integral-derivative (PID) temperature con- 
trollers. The temperature of each sapphire window can 
be controlled independently with an absolute accuracy 
of 0.01 K and a relative RMS stability of 0.001 K over 
24 hours. The annular Peltier elements have an inner 
radius of r = 6.5 mm that determines a clear aperture 
suitable to perform optical measurements. The perfor- 
mances and reliability of the cell have been established 
during a long series of experiments on devices sharing 
a similar conceptual design (see [5T1 [25] and references 
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FIG. 1. (color online) Shadowgraph images of the convec- 
tive patterns: Ra = 1320, a = 24mrad a, Reference image 
taken at t = 0s. b-f, Difference images; dark zones denote 
the warmer fluid, which is less concentrated in the denser 
isobutylbenzene component. Different times shown are: b 
150s, c 400s, d 4000s, e 20500s, and f 29000s. The direction 
of the marginal gravitational acceleration parallel to the cell 
plane is indicated by an arrow in b. The directions of motion 
of the lines for the Super-Highway convection are indicated 
by red and blue arrows in e. The size bar in f corresponds to 
the sample thickness h — 1.3mm. 



therein). A h = 1.3 mm thick layer of the mixture is 
hosted inside the thermal gradient cell that can be tilted 
by using a calibrated screw. We estimate that in the 
horizontal configuration the accuracy of the tilt angle is 
better than 2 mrad. The diagnostic method of choice is a 
shadowgraph, a visualization technique widely employed 
in fluid dynamics 23-25]. Its implementation comprises 
a super-luminous diode (Superlum, Broad Lighter S680) 
with a wavelength of A = (683 ± 9) nm, coupled to a 
single mode optical fiber as a light source. The diverging 
beam out of the fiber is collimated by using an achro- 
matic doublet lens of focal length / = 150 mm posi- 
tioned at the focal distance from the lens. No other lens 
is used after the sample cell. A Charged Coupled De- 
vice sensor (Vosskiihler, CCD4000) with a resolution of 



2048x2048 pixels of 7.4 x 7.4 /im 2 is placed at a distance 
of z = (260 ± 10) mm from the detector [22] . 

A typical measurement sequence involves the rapid im- 
position of a temperature difference by heating from be- 
low. As a result, a nearly linear temperature profile is es- 
tablished across the sample in a time tt — h 2 /Dt — 20 
s. During the process, we grab continuously shadow- 
graph images of the convective patterns. The frame rate 
is set to 1Hz during the initial fast kinetics and reduced 
to l/60Hz for the subsequent slower phase. A sampling 
from a typical image sequence corresponding to an incli- 
nation of 24 mrad and to a Rayleigh number Ra = 1320 
(AT = 8.5 K) is shown in Fig. 1 (see also Supplementary 
Movie SI). A reference background image is taken be- 
fore applying the temperature gradient to the fluid mix- 
ture (Fig. la). This image is then subtracted to all 
the images subsequently collected at generic time t. The 
imposition of the temperature difference is followed by 
a latency time where no pattern formation occurs (Fig. 
lb). This featureless phase ends with the appearance of 
convective rolls (Fig. lc). After a few thousands seconds 
the rolls rearrange towards a more ordered quasi-squared 
pattern (Fig. Id), a configuration typical of solutal con- 
vection in binary mixtures with positive Soret coefficient 
heated from below p^HT7] , 

Quite surprisingly, after a time of the order of 10 4 s the 
square patterns start to drift and evolve gradually into 
sequences of columns of liquid moving in opposite direc- 
tions along parallel lanes arranged in a Super-Highway 
configuration (SH) (Fig. le). Eventually, the SH pat- 
terns might fade away and the convective pattern can be- 
come almost featureless (Fig. If). At Rayleigh numbers 
of the order of Ra SH « 1400 (AT = 9 K) the SH pat- 
terns become stable (Supplementary movie S2). At larger 
Rayleigh numbers, the patterns display marked oscilla- 
tions between a drifting square pattern (DS) and an SH 
one. The long time behavior of the patterns is strongly 
affected by the marginal inclination angle and by the 
Rayleigh number (Fig. 2). The phase diagram displays a 
stable Large Scale Flow (LSF) region at small Rayleigh 
numbers, where convective patterns are destroyed by the 
long-term effect of the inclination. At higher Rayleigh 
numbers the columnar convection is more efficient than 
the LSF. Under these conditions stable drifting square 
patterns are present. The intermediate region is charac- 
terized by patterns where the competition between LSF 
and solutal convection becomes apparent and the system 
displays marked oscillations between drifting square pat- 
terns and longitudinal rolls or SH patterns. In a narrow 
region of the phase space (red stars in Fig. 2) the balance 
between LSF and solutal convection gives rise to a stable 
SH convection state in which the columns organize into 
parallel lanes drifting in opposite directions (Fig.le and 
Supplementary Movie S3). The white spots, which rep- 
resent descending columns of colder liquid crossing the 
cell perpendicularly to the plane of the figure, align into 
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FIG. 2. (color online) Phase diagram of long-term convec- 
tive behavior: phase diagram of the convective pattern as a 
function of the Rayleigh number and the inclination angle. 
Symbols represent the different long-term convective states: 
(V) for Large Scale Flow (LSF);(*) for Super-Highway (SH) 
convection; (o) for oscillations; (□) for drifting square pat- 
terns (DS). The two lines mark the approximate boundaries 
of the transition region. 

parallel lines and move in the direction of the residual 
gravitational forces g r parallel to cell surface, while the 
dark spots align into another set of lines and move in the 
opposite direction. The two sets of lines are staggered 
across the cell plane. 

In order to perform a systematic investigation of the 
evolution of the convective patterns we determined from 
each image sequence the time evolution of the contrast 
C(t) of the images. The image contrast is defined as 
C(t) =< [i(t) — i(t) } 2 >s where < ... >g represents 
the average over the pixels of an image, i(t) = I{t)/ < 
I(t) >s is an image normalized by its spatial average, 
and i(t ) is a normalized background image at time t 
before applying the temperature gradient. The contrast 
provides a quantitative estimate of the strength of the 
temperature and concentration modulations generated 
by the convective motions. Under all the explored exper- 
imental conditions, after the imposition of the tempera- 
ture difference the contrast remains constant for a latency 
time t* lasting from tens to hundreds of seconds (Fig. 3) . 
The latency period is then followed by a rapid growth of 
the contrast, and by relaxation oscillations leading to a 
value that remains fairly stable for a second latency time 
tLSF related to the inclination of the cell. 

Interestingly, up to the time tLSF the evolution of 
the contrast under different experimental conditions is 
qualitatively similar. After that time patterns becomes 
strongly affected by the Large Scale Flow that determines 
a dramatic differentiation of the contrast. At Ra = 1320 
the drop of the contrast marks the appearance of almost 
featureless patterns associated to the LSF (black line in 
Fig. 3). The fact that the contrast does not drop back 
to its original value is the signature of the presence of 



barely detectable fluctuations induced by the shear mo- 
tion. At Ra — 1400 the balance between LSF and solu- 
tal convection gives rise to a SH convective regime (red 
line in Fig. 3) with a stable contrast, accompanied by 
small number fluctuations determined by the entrance 
and exit of the columnar structures into the field of view. 
At Ra — 1480 the competition between LSF and solutal 
convection gives rise to the transition between different 
patterns and the contrast exhibits pronounced oscilla- 
tions (blue line in Fig. 3). Finally, at Ra = 1550 solutal 
convection gives rise to square patterns, but the influence 
of the LSF still determines a drifting of the patterns that 
gives rise to small number fluctuations (magenta line in 
Fig. 3). 

The presence of relaxation oscillations is a signature 
that the onset of convection is determined by the desta- 
bilization of concentration boundary layers (BLs). Ba- 
sically, the BLs grow uniformly by diffusion until they 
reach a critical thickness S* beyond which convection sets 
in. An estimate of the critical thickness of the boundary 
layers can be obtained from the critical dimensionless la- 
tency time t* = t*D/h 2 |2l)H2"8"] . where D is the diffusion 
coefficient. Figure 4a-b shows the dimensionless latency 
time r* and critical BL thickness S* , relative to h, plotted 
as a function of Ra. Interestingly, the data in Fig. 4a-b 
fall onto the same curve independently of the inclination 
angle. This confirms that the initial stages of pattern 
formation are not significantly influenced by a marginal 
inclination of the sample. 

After a time tlsf of the order of 10 4 s (Fig. 3) the pat- 
terns start to drift, thus marking the start of the influence 
of the large scale shear flow induced by the inclination 
of the cell. This time is strongly affected by the inclina- 
tion angle, as it is apparent form the lack of overlap of 
the curves for tlsf(Rci) corresponding to different inch- 
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FIG. 3. (color online) Image contrast: contrast of shadow- 
graph image sequences C(t) plotted as a function of time. The 
layer of liquid is tilted at an angle of 24mrad. The curves cor- 
respond to different imposed temperature differences. Contin- 
uous lines stand for Rayleigh numbers: (from top to bottom) 
fla=1550, 1480, 1400, 1320. 
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FIG. 4. (color online) Convection parameters: a, Dimension- 
less critical time r* for the onset of convection, b, Criti- 
cal thickness S* of the dimensionless concentration boundary 
layer, c, Time for the manifestation of the large-scale shear 
flow (LSF) related to the inclination, tlsf- d, LSF time mul- 
tiplied by the inclination angle, t^sfc, the blue line depicts 
the 1/Ra behaviour. Different symbols stand for different in- 
clination angles: (□) are for 6 mrad, (A) for 18 mrad, (o) for 
24 mrad and (V) for 36 mrad. 



as shown in Fig. 4d. Interestingly, the data does not fol- 
low the contemplated (Ra) -1 scaling. For large Rayleigh 
numbers close to the threshold Ra c the product shows 
signs of a divergence. At such large Ra a slowing down 
of the shear flow occurs, indicating the dominance of the 
solutal columnar convective mass transfer over the shear 
flow. 

In the present work we limited our analysis to one liq- 
uid mixture with a positive Soret coefficient; it remains 
an open question whether a similar mechanism can be 
generally observed in other systems where a temperature 
and composition stratification coexist. The most inter- 
esting case is probably the thermohaline circulation in 
oceans [101 HI]. Here the cooling determined by wind 
at the surface of the ocean determines a decrease of the 
temperature and an increase of the salinity through evap- 
oration, in the presence of a variable landscape at the 
bottom of the ocean. The combined effect of these fac- 
tors leads to the large-scale thermohaline circulation, a 
phenomenon that is still not understood well, due to the 
simultaneous presence of other effects that contribute to 
the oceanic currents. Our experimental model system al- 
lows to isolate and study independently the fundamental 
mechanism driving the thermohaline circulation process. 
This will facilitate progress towards a fundamental un- 
derstanding of the interplay between large scale flow and 
local columnar flows. 
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nation angles (Fig. 4c). A simple understanding of the 
mechanism behind the influence of the shear flow can be 
attained by neglecting the effect of convection (2i. The 
inclination provides a component of the buoyancy force 
Ffc cx gApTSma(l + ip) aligned parallel to the bound- 
aries, where a is the inclination angle and ip = Ap c /Apx 
is the so-called separation ratio. This force determines 
a laminar flow of the boundary layers in opposite direc- 
tions that drives the large-scale flow inside the cell. The 
buoyancy force is contrasted by the shear force Fs oc r]is, 
where v oc 1/t l sf is the mean velocity of the drifting 
boundary layers. At steady state the two forces balance 
each other and, as Ra oc Ap, for small inclination angles 
t lsf oc (Ra) -1 . This simple relation suggests that the 
product tlsfo, all data should collapse on a single mas- 
ter curve, which is confirmed by our experimental results 
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